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Abstract- Differential equations are encountered very often in 
engineering problems and generally in al sciences. Modeling the 
effect of variation of physical quantities such as temperature, 
pressure, velocity, stress, strain, current moisture and many 
other on engineering problems requires most of the times the 
establishment of differential equations. For simplicity reason 
the parameters and variables involved which are measured or 
estimated from experience are considered exact even though 
they often contain uncertainties. One way do deal with these 
uncertainties nowadays is through convex fuzzy sets. According 
to all the above it is almost unavoidable to introduce fuzzy 
parameters and variables in the solution of differential 
equations. Much research was carried out during the recent 
years in theoretic and applied subjects containing fuzzy 
differential equations with H-derivative. This method though, 
in some cases has some disadvantages leading to solutions with 
increasing support as time t increases. In order to alleviate this 
disadvantage the generalized differentiability (G-H derivative) 
was introduced. In this paper the case of a semi-confined 
aquifer is studied, which is bounded on top by a thin 
semi-permeable layer and on bottom by an impermeable layer. 
This system leads to a second order differential equation with 
fuzzy boundary. The solution of this problem is obtained using 
the generalized H-derivative. 

Index Terms — Fuzzy applications, fuzzy differential 
equations, semi-confined aquifer. 


I. Introduction 

Differential equations are encountered very often in 
engineering problems and generally in al sciences. Modeling 
the effect of variation of physical quantities such as 
temperature, pressure, velocity, stress, strain, current 
moisture and many other on engineering problems requires 
most of the times the establishment of differential equations. 
Also, the impact of certain physical quantities on other 
physical quantities leads to differential equations. So, the 
theory of differential equations is found in mechanical 
vibration or structural dynamics, heat transfer, hydraulics etc. 
The engineers should be able to model actual problems using 
mathematical equations and through their solution to 
comprehend the behavior of the systems under consideration. 
Differential equations are extensively used in mathematical 
modeling and engineering applications. Theory and 
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techniques for solving differential equations are applied to 
solve practical engineering problems. 

For the sake of simplicity the parameters and variables 
involved in the systems are considered as crisp or defined 
exactly. Real problems though lack exact information about 
the variables and parameters obtained by experiment or 
experience. Nowadays, these uncertainties are modeled 
through convex normalized fuzzy sets. So, it is necessary to 
include in the solution of differential equations fuzzy 
variables and fuzzy parameters. 

Fuzzy set theory is a powerful tool for modeling uncertainty 
and for processing subjective information in mathematical 
models. The development of this subject has been focused to 
theoretical topics ([ l]-[5]-[14]-[16]) as well as to applications 
like population models, civil engineering and hydraulics 

([11H12]). 

Initially differentiable fuzzy functions were studied by [18], 
who generalized and extended the concept of [13], 
differentiability (H-derivative) of set valued mappings to the 
class of fuzzy mappings. Also, [14] and [19] developed a 
theory for fuzzy differential equations. 

Many works have been carried out during last years in 
theoretical and applied topics for fuzzy differential equations 
with H-derivative ([14]-[16]-[19]-[21]). But in some cases 
this method suffers certain disadvantages that lead to 
solutions with increasing support as time t increases ([5]-[9]). 
This proves that in some cases this solution is not a good 
generalization of the corresponding crisp case. In order to 
surpass the above deficiency, the generalized differentiability 
was introduced ([2] -[3] -[4). This new derivative is defined 
for a larger class of fuzzy functions than Hukuhara derivative 
[13]. 

In this paper the case of a two-point fuzzy boundary value 
problem for a second order fuzzy differential equation is 
studied, the problem concerns the case of a semi-confined 
aquifer which is bounded on top by a thin semi -permeable 
layer and on bottom by an impermeable layer. This system 
leads to a second order differential equation with fuzzy 
boundaries. The solution of this problem is obtained using the 
generalized H-derivative. 

H. MATHEMATICAL MODEL 

A. Definitions 

Definition 1. A fuzzy set U on a universe set X is a 
mapping U ! X — » [0,1] , assigning to each element XeX 
a degree of membership 0<U(x)<l. The membership 
function is also defined as ftp (x) with the properties: 
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|i^ is upper semi continuous, 

ft ~ (x) =0, outside of some interval [c, d], 

there are real numbers a and b, C<&<b <d such that 


is increasing on [c,a], decreasing on [b, d] and [i ~ (x) =1 
for each x G [a,b]. 

U is a convex fuzzy set (i.e. 

H 0 (Ax + (1 - A)x) > min{ fi 0 (Ax), n 0 ((1 - A)x)} 


Definition 2. Let X being a Banach space and U being a 
fuzzy set on X. We define the a-cuts of U as 


[U] a = {x G X U(x) > 0} , with 0 < a < 1 , and for a=0, 


we also define the closure as 


[U]° = {xeX U(x) > 0} . 


Definition 3. Let K(X) the family of all nonempty compact 
convex subsets of a Banach space. A fuzzy set U on X is 
called compact if [U] a G K(X), Va G [0,1]. The space of 
all compact and convex fuzzy sets on X is denoted as F (X). 
Definition 4. Let UgF (R). The a-cuts of U , are: 

[U]“ =[u;(x), u;(x)]. 

According to the representation theorem of [15] and the 
theorem of [10] the membership function and the a-cut form 

of a fuzzy number U are equivalent and in particular the 
a-cuts [U] a = [U~(x), U+(x)] uniquely represent U , 
provided that the two functions are monotonic 
(U a increasing, U* decreasing) and Uj < U^, for 
a=l. The derivatives of each a-cut with respect to x, for a 
given a G [0,1], are denoted as: 


■ _ au;(x) . _au;<x) 

l U al x a ’ l U al x “ 3 

OX OX 

Definition 5. The following notations are used for the above 
derivatives with respect to a 

s (( u;)\) = A(£W) = i.(«W) = 

da dx ox d a 

= (SU“) x(x), 

5((u:yo= A ( Mw )= A ( Mw )= 

da ox ox da 

= (su;) , x(x) 

Lemma 1. If both U a (x), and U^(x) are differentiable 

w.r.t. x, then the a-cuts of the gH (generalized 
Hukuhara)-derivative of U are: 

U gH (x) = [mm{(U a )'(x),(U;)'(x)}, 

max{(U”) (x), (U+) (x)}] 

provided that the two functions: 

Ug H (x)“ =min{(U“) (x), (U a + ) (x)}, ^ 
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U gH (x); =max{(U”) (x), (U;)'(x)}, 
define a fuzzy number (w.r.t. a). Here the functions 
(U a ) , (LT) denote the derivatives with respect to x, for 
given a G [0,1] . 

Definition 6. Let [U] G F(R), with the 

a-cuts [U] a = [U~(x), U^(x)] satisfying Lemma 1. 
According to [3] : 

1. U is (i)-generalized Hukuhara (gH) differentiable if: 

(u-)x(x)<(un'x(x) 

(6u;),(x) = |-(|-u;(x))>o 

da dx 

(8u;)’,(x) = -f c-f u;(x))so 

da dx 

for Va e [0,1]. 

2. U is (ii)-generalized Hukuhara (gH) differentiable if: 

(U-)'x(x)>(U0'x(x) 

(6U-)x(x) = A ( A u - (x)) <o 

da dx 

(su;Ux) = j-(j-u;(x»zo 

da dx 

for Va g [0,1]. 


B. Seni-confined aquifer 

The semi-confined aquifer (Fig.l) is discharging towards a 
Lake. This aquifer is bounded on top by a thin 
semi-permeable layer of thickness B ' and on the bottom by an 
impermeable layer. On top of the semi-permeable layer there 
is water entrapped by an earth dam. The equation describing 
the water level variation h as a function of length x is given 
as: 


d 2 h h -h 

dx 2 + A 2 


= 0, 


(1) 


where A, is a leakage factor given by the relation: 


A = a/KBc , (2) 


and c is the hydraulic resistance of the semi -permeable 
aquifer having time dimensions and is given by the following 
relation: 


C = 


K 


(3) 


Equation (1) resulted from the application of Darcy’s law and 
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also taking into account the leakage from the 
semi-impermeable layer ([19]-[21]). In relations (1), (2), (3), 
K is the hydraulic conductivity of the semi-confined aquifer, 
B is its thickness, B' is the thickness of the thin 
semi-permeable layer , K' is its hydraulic conductivity, hO is 
the head of the aquifer for X — » oo , and hi is the head of the 
aquifer at x=0. 



s(X,a) = Cje x/x +c 2 e 


-x/A, 


For X=0, s 0 =Ci+c 2 , and for X — » oo, C, = 0. 
So c 2 = s 0 and the solution now takes the form: 


s(X,c) = s 0 e x 


h = h 0 - s 0 e“ 


or, 


D. Fuzzy model 
We fuzzify now the solution: 
h = h 0 - s 0 e _x , (10) 


(9) 


Fig. 1. Semi-confined aquifer. 

We introduce now in (1) the drawdown s=ho-h and we 
obtain: 


dx 2 (4) 

The boundary conditions of the above problem then become: 

s(0) =Sj =h 0 -h p x = 0 

s(oo) = s G = h 0 -h 0 = 0 x^co. 


We replace now in (4): 



A 

5 


and this leads to : 



with boundary conditions : 

s(0) =s 0 =h 0 -h p X = 0 

s(oo) = s m = h 0 - h 0 = 0 X— »oo. 


( 6 ) 


(7) 


( 8 ) 


C. General solution 

The above equation has the general solution: 


where s 0 = (r, ir^), = (0, m 2 ), and ml, m2 are the 

centers of fuzzy numbers and r, 0 are the spreads. The number 
S 0 is a symmetrical triangular fuzzy, while , ( m 2 =0), 

is a crisp number, expressed in fuzzy notation. The fuzzy 
boundary conditions are now: 

s 0 |“ =[s a <rf,s a or] =[ 111 ! -r(l-a), m! +r(l-a)], 

Soo | “ = [s a cor , s a oor ] = [0, 0] , (m 2 = 0) 

From (9) and (10) we arrive to: 

h|“ = h 0 | a -^| a e" x =[h 0 ,h 0 ]-[Soe' x , Sge" x ] = 

= [(h 0 -s;e- x ),(h 0 -s 0 e^ x )] = 

=[{h 0 -(m!+r(l-a))e~ x }, 

[h 0 - (m! - r(l - a))e“ x } = [ h(a, X) - , h(a,X)] . (J2) 

The functions h(a, X) and h(a,X) + are both 
differentiable w.r.t. x and we have: 

h x (a,X)“ = (m 1 +r(l-a))e“ x , h' x (a,X) + = 

=( m ! -r(l-a))e” x ? (13) 

where 

(mj + r(l - a))e x > (rtij - r(l - a))e x , or 
h x (a, X)“ > h x (a, X) + . 

According to Lemma 1 the gH-derivative of h is: 

[h]' x | a ={min[h' x (a,X)', h' x (a,X) + ], 

max[h x(a,X) _ , h x (a,X) + ]} = 

= [min{(m! +r(l-a))e" x , (m, -r(l-a))e“ x }, 
max{(m! +r(l-a))e“ x , (m, - r(l-a))e" x }] = 
= {h (x)]“,[h (x)]^}, 
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providing that the two functions: 

h (x)|” =mm[h x (a,X)“, h x (a,X) + ] 

1 and 

h( x )|a =max[h x(a,X)“, h x (a,X) + ] 

define a fuzzy number with respect to a. 

We have now: 


Fig. 2. Boundary conditions 


E. Discharge calculation 


The discharge per unit length is given by: 

dx A dX (i 5 ) 


h (x)| a =min[h x (a,X) , h x (a,X) + ] = 

= m, - r(l - a)e x , 

h(x)| + = max[h x (a,X)“, h' x (a,X) + ] = 

= nij + r(l - a)e” x 

5h' x (a,Xr _ lm _x 5h x (a,X) + _ x 

— “I - re y — re 

da da 

h'(x)|"=h'(x)|^ 

So the two functions h (x)|~, h (x)| + , define a fuzzy 
number with respect to a. 

For the above relations we have: 

h’ x (l,X)-=m 1 =h x (l,X) + = m 1 

(6h;)x(x) = ^-(^-f;<x)) = 
da ox 

_ 0{mj + r(l - a)e“ x } _ 

da 

= -re“ x < 0, for Va e [0,l], (14a) 

(8h;)' x (x) = ^-(^-f a + (x)) = 
da dx 

_ dfrrij - r(l - a)e~ x } _ 

da 

= +re~ x >0, for Vae [0,1] , (J4b) 

According to [3], (14a) and (14b) determine that the function 
h(a, X) is (ii) gH-differentiable. 


The derivative dh/dX is taken from (13) and (15) takes 
now the following form in fuzzy notation: 


q “ = — [(m, - r(l - a))e x , (m, + r(l - a))e^ x ] 

A 


HI. APPLICAIONS 


A. Application 1 

It is assumed that the aquifer has the following 
characteristics: 

K=0.4cm/s (gravel), K'=0.01cm/s (sand), hl=B=20m, 
h0=24m, hl=20m, B"=0.8m, X=V( KBB7K')= 25.29m, 
KB/X=0.003162m/s, ml=4m, r=lm. 

The solution is: 

h|“ =[{h 0 - (nij +r(l-a))e“ x ], 


24 

23 

22 


21 
20 
19 

0 1 2 3 4 5 6 

Fig. 3 Variation of solution h(a,X) as a function of 
distance X. 


[h 0 - (nij -r(l-a))e” x }] = 
= [ [24 - (4 + (1 - a))e” x } , 

{ 24 - (4 - (1 - a))e” x } ] 


h(a,X 
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Fig.4 Membership function of the solution at x/k=0, 1,2. 
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Fig.5 Variation of q as a function of X. 
In position (X=0) the discharge towards the Lake is: 


h|“ =[{h 0 -(nij + r(l-a))e” x }, 
{h 0 - (nij - r(l - a))e” x } ] = 
= [ {22 - (6 + (1 - a))e” x } , 

{ 22 - (6 - (1 - a))e“ x } ] 
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° ] =^ig-7. Variation of h(a, X) as a function of distance X. 


= 0.003 162 [(a + 3), (5 - a)] 

This discharge for a=0 , takes the following values: 

q":° 0 =[0.0095, 0.0158 JmVnT 1 . 

Assuming that the dam has a crest lengthlOO m, the discharge 
towards the Lake takes the following value: 

q“:° 0 = [0.948, 1. 581 ]mV, 

the discharge volume will be: 

V“2 = [81966, 136 610]m 3 /d, 

the volume for a=0.5 takes the following value : 

V“:|f = [95 627, 122 949]m 3 /d. 

and finally for a=l the discharge volume will be: 

V“1 = [109 288, 109 288]m 3 /d. 


’ and for time span of one day 
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Fig.6 Discharge volume V towards the Lake. 

B. Application 2 

We assume an aquifer with the following characteristics: 
K=0.2cm/s(gravel), K' = 0.04 cm/s(sand), B=16m, 
B' = 3 m, ^=1 5.49 m, ml=6m, l-l m, s(0)=6m, 

KB/X=0.002066m/s, h0=22m, hl=16m. 
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In position (X=0), the discharge towards the Lake is: 

k 

(m 1 +r(l-a))e“ x=0 ] = 

= 0.002066 [(a + 5), (7 -a)]. 

This discharge for a=0, takes the following value: 

q“:° 0 = [0.01033 , 0.01446 ]m 3 s" 1 m" 1 . 

Assuming again a dam with 100m crest length: 

q“:° 0 =[1.033, 1.446 ]m 3 s _1 , 

and during a day the discharge will be: 

V“2 =[89251, 12495 l]m 3 /d. 

the volume for a=0.5 takes the following values: 

V“2 5 =[98157, 116003]m 3 /d. 

and finally for a=l the discharge becomes: 

V“^ =[107101, 107101]m 3 /d. 
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Fig-8 Variation of q as a function of X. 


The solution is: 
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Fig.9 Discharge volume V towards the Lake. 

C. Application 3 

We assume an aquifer with the following characteristics: 
K=0.2cm/s(gravel), K'=0.01cm/s(sand), hl=B=20m, 

h0=24m, B"=0.8m, X=V( KBB/K> 17.88m, 

KB/AM). 002236m/s, ml=4m, r=lm. 

The characteristics are the same as in the example 1, except 
the hydraulic conductivity K, which is the half of the value of 
example 1. For that reason the values of X and KB/X have 

changed. The variation of h(a, X) vs. X, remains the same 
as in Fig. 1, with different X, but the variation of q vs. X, as 

well as the value of discharge volume V towards the Lake 
have completely changed, because of the change of X and the 
KB IX. The fig. 10, 11, 12 show the existing differences 
between the example 1 and example 3. 
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Fig.10 Variation of h(a, X) as a function of distance X. 
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Fig. 12 Discharge volume V towards the Lake. 

IV. COMMENTS-CONCLUSION 

The function h(a,X) is (ii) gH-differentiable and its 
derivative h x (a, X) equals: 

[h] x | a = {mn[h x (a,X)“, h x (a,X) + ], 

max[ h x (a,X)“, h x (a,X) + ]} 

and so a solution to this problem can be obtained. 

The fuzziness is decreased with increasing X=x/X and 
practically for X > 4 in both cases the values of h and q 
became crisp numbers. 

As it is shown in the above examples the discharge volume 
decreases as a, which is called the confidence level, is 

increased. Consequently the interval [ V r a , V^ a ] is decreased 

by 50% accordingly to the increase of a. 

The discharge volume in the second case is decreased by a 
percentage of 2% compared to the volume in the first, even 
though the hydraulic conductivity is reduced by 50% and the 
thickness of the aquifer is reduced by 20%. This is due to the 
fact that the value of KB/A, in second case multiplied by the 
slope of the head at x=0, provides a value which is only 2% 
less than in the first case. 

The discharge volume in the third case is decreased by a 
percentage of 30% compared to the volume in the first, even 
though the hydraulic conductivity is reduced by 50%. This is 
due to the fact that the value of leakage factor X in third case is 
reduced by 30%, as well as the value of KB/X , and q is 
multiplied by this value. 

So in our case the leakage factor and the slope of the head 
function have a very significant role in the solution. 
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Fig.ll Variation of q as a function of X. 


REFERENCES 

[1] Abbasbandy S., J.J. Nieto, M. Alavi, “Tuning of reachable set in one 
dimensional fuzzy differential inclusions”, Chaos, Solutions and 
Fractals, vol. 26, 2005, pp. 1337-1341. 

[2] Bede B. and Gal S.G., “Generalizations of the differentiability of 
fuzzy-number- valued functions with applications to fuzzy differential 
equations”, Fuzzy Sets and Systems, vol. 151, 2005, pp. 581-599. 

[3] Bede B. and Stefanini L., “Solution of Fuzzy Differential Equations 
with generalized differentiability using LU-parametric representation”, 
EUSFLAT-LFA , 2011, pp. 785-790. 

[4] Bede B. , “A note on two-point boundary value problems associated 
with non-linear fuzzy differential equations”, Fuzzy Sets and Systems, 
vol. 157, 2006, pp. 986-989. 




le\t$en 

Research Publication 


28 


www.ijntr.org 


gy 






































4 

























l 

59251 / 10‘ 

noi 

1245 

>51 

V m 3 /d 


q(m 3 /s 

/m) 























\ q, 

tx.oi 







\V 





A=17.8 

8m 



\ c 

b(X,l) 

= q r (X, 

1) 





A' 

\ 







q,(X,0 

) ^ 







X=x/A 




7 

r 



‘ V 
















\ 











* 





T 
















57959/ 

77276 


96595/ 

V m 3 /d 


International Journal of New Technology and Research (IJNTR) 
ISSN:2454-4116, Volume-2 , Issue-9 , September 2016 Pages 23-29 


[5] Bhaskar T., V. Lakshmikantham, V. Devi, “ Revisiting fuzzy 
differential equations”, Nonlinear analysis, vol. 58, 2004, pp. 351-358. 

[6] Bucley J .J and T. Feuring, “Fuzzy differential equations”, Fuzzy Sets 
and Systems, vol. 110, 2000, pp. 43-54. 

[7] Bucley J.J and T. Feuring, “Fuzzy initial value problem for Nth-order 
linear differential equations”, Fuzzy Sets and Systems, vol. 121, 2001, 
pp. 247-255. 

[8] De Ridder N. A. and G. Zijlstra, “Seepage and Groundwater Flow”, in 
Drainage Principles and Applications, ILRI Publication 16, H.P. 
Ritzema Ed., 1994, pp. 305-340. 

[9] Diamond P., “Brief note on the variation of constants formula for fuzzy 
differential equations”, Fuzzy Sets and Systems, 129, 65-71, 2002. 

[10] Goetschel R. and W. Voxman, “Elementary fuzzy calculus”, Fuzzy 
Sets and Systems, vol. 18, 1986, pp. 31-43. 

[11] Guo M., X. Xue, and R. Li, “The oscillation of delay differential 
inclusion and fuzzy biodynamics models”, Mathem. Comput. 
Modelling, vol. 37, 2003, pp. 651-658. 

[12] Guo M., X. Xue, and R. Li, “Impulsive functional differential inclusion 
and fuzzy population models”, Fuzzy Sets and Systems, vol. 138, 2003, 
pp. 601-615. 

[13] Hukuhara M., “Integration des Apllications Mesurables dont la Valeur 
est un Copact Convexe”, Funkcialaj Ekvacioj, 10, 205-233, 1967. 

[14] Kaleva O., “Fuzzy differential equations”, Fuzzy Sets and Systems, 
vol.. 24, 1987, pp. 301-317. 

[15] Negoita C.V. and D.A. Ralescu, “Representation theorems for fuzzy 
concepts”, Kybernetes, vol. 4, Iss 3, 1975, pp. 169 - 174. 

[16] Nieto J.J. and R. Rodrfguez-Lopez, “Bounded solutions for fuzzy 
differential and integral equations”, Chaos, Solutions and Fractals, vol. 
27, 2006, pp. 1376-1386. 

[17] O’Regan D., V. Lakshmikantham, J. J. Nieto, 2003, “Initial and 
boundary value problem for fuzzy differential equations” Nonlinear 
analysis, vol. 54, 2003, pp. 405-415. 

[18] Puri M.L. and Ralescu D. A., “Differentials of Fuzzy Functions”, 
Journal of Mathematical Analysis and Applications, vol. 91, 1983, pp. 
552-558. 

[19] Seikkala S., “On the fuzzy initial problem”, Fuzzy Sets and Systems, 
vol. 24, 1987, pp. 319-330. 

[20] Terzides G. and D. Karamouzis, “Ground Water Hydraulic”, Eds Ziti, 
Thesaloniki-Greece, 1985, pp. 75-80. 

[21] Vorobiev D. and S. Seikkala, “Towards the theory of fuzzy differential 
equations”, Fuzzy Sets and Systems, vol. 125, 2002, pp. 231-237. 



Christos Tzimopoulos 

Employment Address: School of Surveying & Rural Engineering , 
Aristotle University of Thessaloniki 54124- GREECE, Tel.: 
+302310996141- Fax:+ 302310996049- E-mail: ctzimop@gmail.com. 
Academic degrees: 1962. M.Sc. (Civil Engineering) Aristotle University of 
Thessaloniki. Thessaloniki-Greece, 1970. Ph.D. (Civil Engineering) Aristotle 
University of Thessaloniki. Thessaloniki-Greece, 1973. Dr.Ing. Universite 
de Paris VI — Paris-France, 1975. Habilitation (School of Agriculture) 
Aristotle University of Thessaloniki, 1975. Assistant Prof. Aristotle 
University of Thessaloniki, School of Agriculture. Thessaloniki-Greece, 
1977-2006. Professor, University of Thessaloniki, School of Rural and 
Surveying Eng. Thessaloniki-Greece, 2006. Emeritus Professor, Aristotle 
University of Thessaloniki. Thessaloniki-Greece, 2014 Member of Greek 
Agricultural Academy. (GAA). 2016 Honorary Doctor of the Democritus 
University (Department of Civil Engineering) of Thrace-Greece. (1985 
Professeur Associe-Ecole Nationale d’Hydraulique de Grenoble, 
I.N.P.-France.,1987 Professeur Invite - Ecole Nationale d’Hydraulique de 
Grenoble, I.N.P. -France. Professional Activities: 1966-1977: During this 
period he has worked as engineer in the Department of Water Resources and 
Hydraulic Works, in the Ministry of Public Works, Greece. Research 
Activity: Hydrology (surface and subsurface), water resources management 


and processing of digital hydrologic conditions, fuzzy engineering 
applications. Publications: He is the writer and co-writer of over of 270 
scientific papers related with water resources, published in Greek and 
International Scientific Journals as well as in International Congresses. He 
was the main supervisor in more than 18 Ph.D theses concerning water 
resources. 



Chris Evangelides 

Employment Address: School of Surveying & Rural Engineering , 
Aristotle University of Thessaloniki 54124- GREECE, Tel.: 
+302310996147- Fax:+ 302310996147- E-mail: evan@vergina.eng.auth.gr. 
Academic degrees: B.Sc in Electrical engineering from California State 
University Long Beach, M.Sc in Electrical Engineering from California State 
University Long Beach, Ph.D in Aristotle University of Thessaloniki School 
of Surveying and Rural Engineering, Associate Professor at the Department 
of Hydraulic Works and Transportation Engineering. Research Activity: 
vadoze zone, irrigation, water resources management, application of fuzzy 
logic in water resources management, optimization, experimental processes 
in infiltration and absorption of soils,. Page in Research Gate : 
www.researchgate.net/profile/Chris_Evangelides 



Kyriakos Papadopoulos 

Kyriakos Papadopoulos obtained his B.Sc. in mathematics from the 
University of Exeter (UK) and his master degrees from 
Liverpool/UK(Mathematical Sciences) and Patras/Greece (Pure 
Mathematics), respectively. He did his Ph.D. in the University of 
Birmingham(UK); his Ph.D. Thesis title is "Set Theoretic and Topological 
Characterizations of Ordered Sets". He is currently an Assistant Professor of 
Mathematics, in the American University of the Middle East, in Kuwait. 



Basil Papadopoulos 

Basil Papadopoulos is a Professor of Mathematics and Statistics at 
Democritus University of Thrace, Department of Civil Engineering, Section 
of Mathematics and Informatics. His recent research topic is Fuzzy Logic 
with Engineering Applications. He has also worked in the past in the 
following research areas: topology and analysis. He has published about a 
hundred papers in research Journals and Conferences. See Google Scholar: 
(https://scholar.google.com/citations ?user= ZlaVyMcAAAAJ&hl= 
en&oi=ao) and personal website: http://utopia.duth.gr/papadob/ 




le\t$en 

Research Publication 


29 


www.ijntr.org 



